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NOMENCLATURE 


Eq speed of sound in undisturbed fluid 
A,B integration constants 

dj^j deformation rate tensor 

h parallel channel height 

I^,l 2 dissipation integrals (see eq. 22)) 

k coefficient of heat conduction 

K vortex strength 

L length of vortex 

M Mach number 

p pressure 

r radial distance 

r^ viscous core radius 

S entropy 

t time 

T temperature 

u,u^ velocity vectors 

Ug tangential velocity component 

u streamwise velocity 

Uq peak jet velocity at the nozzle exit 
W parallel channel width 

,x^ position vector 

X streamwise direction 

y direction normal to parallel channel wall 

tnaeam net bunk nn 


+ 3 


Kronecker delta 


^iJ 

X second coefficient of viscosity 

n 3Ug/3r 

E, 3u/3y 

()) velocity potential 

$ viscous dissipation function 

radiative component of viscous dissipation function 

incompressible approximation to dissipation function 
p density 

V kinematic viscosity 

y viscosity 

vorticity 

( ) time average of the quantity in parentheses 

( )• time-dependent fluctuation of the quantity in parentheses 
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VORTICAL DISSIPATION IN TWO-DIMENSIONAL SHEAR FLOWS 


W. Clifton Horne and Krishnamurty Karamcheti* 
Ames Research Center 


SUMMARY 


An exact expression is derived for the viscous dissipation function of a real 
homogeneous and isotropic fluid, which has terms associated with the square of 
vorticity, wave radiation, and dilatation. The implications of the principle of 
maximal dissipation rate, as expressed by Ziegler and others, are explored by means 
of this equation for a parallel channel flow and a cylindrical vortex flow. The 
consequences of a condition of maximum dissipation rate on the growth of distur- 
bances in an unsteady, laminar shear layer are apparently consistent with predic- 
tions and observations of maximum growth rate of vortical disturbances. Finally, 
estimates of the magnitudes of several dissipative components of an unsteady vortex 
flow are obtained from measurements of a periodic wall jet. 


INTRODUCTION 


The exact conservation equations of mass, momentum, and energy, which govern 
the motion of a real fluid are complete, yet intractable for most flows of interest, 
particularly for cases of unsteady or turbulent flow. Attempts to obtain solutions 
for the time-averaged form of these equations involve additional difficulties in 
closure because of the average product terms. The nonlinear terms in these equa- 
tions have required the imposition of simplifying assumptions and approximate 
models, depending upon the particular flow under investigation. 

In addition to the conservation equations, the second law of thermodynamics, 
which requires a net increase in entropy for a real process, has been applied to the 
analysis of shock waves and other irreversible phenomena. The phenomenon of turbu- 
lence is apparently consistent with the second law, as many commonly observed flows 
exhibit states of increasing disorder in the streamwise direction. However, a wide 
variety of flows are characterized by the persistence of a global structure such as 
waves and vortex arrays in the absence of global boundary conditions. The organi- 
zation of such flows, many examples of which are presented by Van Dyke (ref. 1), do 
not exhibit simple structural organization characterized by increasing disorder in 
the streamwise direction. 
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Ziegler (refs. 2 and 3) and others (refs. 4-5) have proposed an extension of 
the second law, termed the principle of maximal rate of entropy production. This 
principle is presented alternatively as the principle of maximal dissipation rate, 
the principle of least dissipative stress, or a principle of least deformation 
rate. The principle states that for a small material "element" throughout which 
properties are uniform, the observed motion results in a maximum of the dissipation 
rate, i.e., the rate of entropy production. This is not to be confused with the 
maximum entropy method of information theory, in which the entropy is associated 
with the loss or lack of information and the maximization of this quantity leads to 
an optimal representation of the power spectrum of a finite sample set (ref. 6). 
Ziegler shows these principles to be the consequence of a thermodynamic orthogonal- 
ity condition relating the dissipative stresses and strains. In a later paper, 
Ziegler (ref. 7) employs the maximum entropy principle to derive constitutive rela- 
tions for heat conduction and deformation in linear elastic and thermoplastic 
solids, and in inviscid and viscous gases and liquids.^ 

Zeigler derives the orthogonality principle from a tensor expansion of the dis- 
sipation function. The associated extremum principles are based on the conditions 
that, for a real substance, surfaces of constant dissipation rate in deformation 
rate space are strongly convex (i.e., a plane which is tangent to a dissipation sur- 
face intersects that surface at one and only one point), and that along any ray 
extending from the origin of deformation rate space, the dissipation rate increases 
in value faster than a linear function of the distance from the origin. Ziegler 
also distinguishes between three types of processes, which are described as elemen- 
tary, compound, or complex, depending respectively on whether the process involves a 
single tensor, multiple tensors which are uncoupled and independent, or multiple 
tensors which are coupled and dependent. The particular case of a heat conducting, 
newtonian fluid is described as a compound process (ref. 2, p. 265), and the orthog- 
onality condition is to be applied separately to the viscous and thermal dissipation 
terms, but not necessarily to their sum, which is the rate of entropy production. 

If it can be determined that the maximum entropy principle applies on a global 
basis, this principle would have a significant role in the determination of overall 
structure in real flows. In this paper, we will derive an exact expression of the 
dissipation function of a homogeneous, isotropic, newtonian fluid which has terms 
associated with the square of vorticity, wave generation, and dilatation. The 
implications of the maximal-dissipation-rate principle are explored by means of this 
equation for a parallel channel flow and a cylindrical vortex flow. Areas in which 
there is consistency between the principle and the notions of the instability of 
laminar shear flows are discussed. Finally, estimates of the magnitudes of several 
dissipative components of an unsteady vortex flow are obtained from measurements of 
a periodic wall jet. 
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Related discussions are presented in references 8-12. 
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DISSIPATION FUNCTIONS FOR NEWTONIAN FLOW 


The energy conservation equation may be expressed in terms of the fluid entropy 


as; 


PT f = . . kV^T 


( 1 ) 


The dissipation function, for a newtonian fluid, ®, is written by Ziegler (ref. 2) 
as : 


. = kdjjdjj . 2„djjd^j 


where 


. 1 a . !!i\ 

■ 2 \»*l »*j/ 


Jeffreys (ref. 13) rewrites the latter expression of equation (2) as: 

3u. 


$ = X 


3Ui ^ r 

35]; * “Is “ 


lj"ij * ^ ax 


i!!il 

J ’■‘ij 


where 


u 


! 3u^ ^'^i\ 

ij ' 


(2) 


(3) 


Hinze (ref. 14) obtains a similar decomposition of the time-averaged dissipa- 
tion function, which is valid only for homogeneous turbulence. The last term in 
(3), (3Uj^/3Xj)(3Uj/axj^), appears in several aeroacoustic wave equations, such as 
that derived by Phillips (ref. 15). Another formulation is obtained by combining 
the terms of the equation: 

- ^ + MV-u = 0 (4) 


where 


„ 3P 3P ’“i 

" = 3t * “i 357 * ^ 357 

1 1 


3u. 3U. , 3 t! 

f , P ou. — 1 . 2E i- 

3t 1 3X^ 3Xj 3X^ 


/3U. 3U A 3U. 

ij ■ ^ 3x. j “^ij^ 3X^ 


(5) 

( 6 ) 
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M = f = 0 for a fluid free of sources of mass or momentum. By expanding the 
terms of equation (4), we obtain: 


. 3u . , / .2 ,2 

i .1 _ If 3 P 3 p 

. 3x. ■ “ p \3x.3x, " ^i^^j 3X.3X, 
J 1 ' J j 1 j 


3u_. 3u, 


2u 




J 3Xj3t 


ifpV 1 

atV “ ”‘l“j 


3u^ 3u . 
3x. 3x, 

1 J 


(7) 


By substituting equation (7) into equation (3), the desired form is obtained: 


- 2y 


« = - (2u * X) 


.2 /3U.\ 


+ (X - 2p) 


3u^ 3u^ 
3x^ 3Xj 


2y I 3^p 3^p „ 3^p 3^p \ p 

p \3 Xj3Xj ■ *^i'^j 3x^3Xj “ 3Xj3t " 2 “ij“ij 


( 8 ) 


This exact expression accounts for three important mechanisms of energy dissi- 
pation in real flows: (1) irreversible expansion or compression, (2) generation and 

radiation of sound or shock waves, and (3) generation of vorticity. According to 
the second law of thermodynamics, the overall entropy should increase for a real 
process, which implies a positive definite value for An extremum condition, as 

proposed by Ziegler on the rate of entropy production is a stronger constraint which 
may guide the determination of the distribution of energy dissipated by the various 
mechanisms . 

The radiative dissipation term: 
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(9) 


may be simplified for comparison with other wave equations by first assuming an 
isentropic relationship between the pressure and density, and then examining the 
wave generation in flows with negligible velocity in all but one direction: 


For these conditions: 



u^ = U >> u^.u^ 


M = 


$ 


R 



- M^)7^p' 


p IL 3^P' _ j_ 3^P' I 

a^ 3x.3t 2 .2 I 

0 1 a 3t I 

o "* 


( 10 ) 


This expression may be compared with the compressible, unsteady, potential flow 
equation: 
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( 11 ) 


M2\n2. 2M 3^(|) 1 3^<1) _ n 

O ci O V 

o 


For waves of acoustic intensity, $ should be negligible except in the sound- 
producing region. For shock waves, the isentropic pressure-density relationship is 
invalid. 


The principle of maximum dissipation rate is established for elemental volumes 
that contain a substance of uniform material properties and velocity. We now seek 
to determine if the principle holds on a global basis for selected flows. As with 
other physical principles or laws, the principle is validated by observations of 
consistency with real flows. 


ENTROPY CONSIDERATIONS IN LOW-SPEED PARALLEL AND 
CYLINDRICAL SHEAR FLOWS 


For incompressible, uniform-density flows, the expression (eq. (8)) for the 
dissipation function reduces to: 


p_ 


7^p 


+ U 0) 


( 12 ) 


According to the principle of maximum dissipation rate, we seek to find the extrema 
of the volume integral of equation (12) for selected familiar flows. For parallel 
channel flow, the mass and momentum equations reduce to: 


2 

3u lE P_ 3 u 

3t ■ ' p 3x p ^ 2 

o o 3y 


For this case, the term in the dissipation function (eq. (12)) associated with 
the Laplacian of the pressure is zero, as is seen by taking the x-derivative of the 
momentum equation and by using the continuity equation: 


Therefore, 



2 2 
i_ 1 

^o 3X^ ^o 3y^ 




3X 


2 


= 0 


(13) 


(14) 
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The volume integral of the dissipation function is now given by: 


I 


dv = y 


/„ 


0 ) dv = uxW 


f'^2 /^Y 

Jy. VV 


dy 


(15) 


where x and W are the length and width, respectively, of the channel. 

To find the extrema of the preceding integral, we seek solutions of the Euler- 
Lagrange differential equation (see ref. 16) 


3$ ^ . . 

3u ■ dy ’ 


5 = - 

3y 


For this example: 


0 . . .2 ^ - 0 

dy2 


( 16 ) 


(17) 


The dissipation extrema are found for the linear profile: 

u(y) = Ay + B 


Thus the linear velocity profile of a parallel channel flow corresponds to a 
dissipation extremum. This profile is observed when the channel walls are in 
relative motion and there is no pressure gradient in the channel. 

Another simple example is that of a cylindrical rotational flow, in which the 
tangential velocity and pressure are functions only of the radius from the origin. 
Hence: 


u = (u ,u„u } ; u = u = 0 
r’ 9 z ’ r z 

The radial pressure gradient is given by the radial momentum equation: 


1 l£ - 

p 3r “ r 


The Laplacian of the pressure is written as: 

„2 


1 d / dp\ 1 d . 2, 

p = F a; (■• = ? d? ‘‘’“e* 


The vorticity is given by: 


, * s - Id, , ““e “e 

“ = y X u = p 5? (rUg) = * r 


( 18 ) 


(19) 
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( 20 ) 


2 2 ^ :;i 

“ ■ Vdr / ^ r dr 2 


We seek the extrema of the integral: 


/• 


f 


dv = - — 2nL I 7^pr dr + y2irL | u^r dr 
P 


a 

/ 


where L is the length of the vortex. The first integral on the right side 
vanishes, since from equation (18): 


J V^pr dr = J p 


d / x2 . 

dF 


r - 


^) dr = 

0 9 


and 


for any realizable flow. Then: 


Ue(o) = Ug(“) = 0 




r dr 


The extrema are obtained as solutions of the Euler-Lagrange equation: 


d_ . 

3u dr \3n/ ^ ’ 


n = 


dr 


* « r 


\dr / ^ ^0 dr ^ r 


du u 

3«_ . 2 i + 2 — 

3u„ ‘dr r 


( 21 ) 


( 22 ) 


(23) 


(24) 
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Then 


3n 


dr 


(s) 


2r 


du 

dr 

du 


0 


2u. 


0 


d^u. 


dr 


+ 2r 


+ 2 


dr 


dr 


By substituting equations (24) and (25) into (23), we obtain: 


2 

r 


d^u. 


dr 




(25) 


This is recognized as Euler's equation, with solutions: 

Ug “ r,r"^ 

These two solutions represent respectively the core and the outer-potential region 
of a Lamb vortex: 

u„ = Kr for 0 < r < r {viscous core} 

0 o 

Kr^ 

u„ = — - for r < r < “ {inviscid outer flow} 

e r o 

The two solutions represent minimum and maximum dissipation conditions in the 
same flow field. Unlike the previous example for the parallel channel flow, there 
are many solutions for the cylindrical vortex, including the Taylor vortex: 

u„ = ^ exp{-r^/4vt} 

® t^ 

The flow configuration obtained by maximum-dissipation considerations, with a 
viscous inner-flow region surrounded by an inviscid outer flow, is consistent with 
the basic notions of boundary layer theory. Further study is needed to determine 
the global validity of the dissipation principle for general laminar flows and for 
steady, high-speed flow. Another aspect of the dissipation principle (to be consid- 
ered later) is the applicability of the principle, on a time-averaged basis, in the 
determination of the overall structure of unsteady and turbulent flow. 


DISSIPATION CONSIDERATIONS IN UNSTEADY, TWO-DIMENSIONAL FLOWS 


The maximum entropy principle is of considerable utility if it can be applied 
as a global constraint for complex flows, and further, if the principle can be 
applied in a time-averaged formulation. Neither of these approaches has been vali- 
dated, and therefore both warrant further investigation. 


8 



The application of extremum conditions of entropy or dissipation rate to 
unsteady flows may be facilitated by considering the mean and fluctuating components 
of the dissipation function: 

p = p + p'; (D = a)-t-(u' 

♦ = - ^ 7^(p + p') + y(« (26) 

1 p 

♦ , = - — ^ 7^p + y[((i))^ + 

1 p 

We may use this expression to interpret the small and large disturbance motions 
in two-dimensional flow. A large number of parallel shear flows such as the 
boundary-layer, the free-shear-layer, and the wall- jet flow have been successfully 
analyzed for stability to small disturbances via the Orr-Sommerfield equation. The 
stability characteristics for these flows are found to be strongly associated with 
the mean velocity profile of the flow. For velocities which exceed a critical 
Reynolds number, small disturbances of the frequency corresponding to the maximum 
amplification rate are predicted to grow exponentially with downstream distance 
until the magnitude of the disturbances exceeds the level permissible in a linear- 
ized analysis. Within the small-disturbance region, the mean velocity profile 
remains unchanged, and the small disturbances take the form of convecting vortical 
motions. Experimental studies of various flows confirm that the frequencies pre- 
dicted for maximum-vortical -disturbance growth rate correspond to the observed 
frequencies of unforced fluctuations. From equation (26), if appears that the 
maximum-vortical-disturbance growth rate corresponds to a maximal dissipation condi- 
tion, since the contribution from the pressure terms is negligible and the contribu- 
tion from the mean vorticity terms is independent of disturbance frequency. 

Downstream of the small-disturbance, exponential-growth region, many parallel 
shear flows exhibit a comparably sized region of periodic large-vortical distur- 
bances. In this region, the mean velocity profile is significantly altered with 
respect to the profile in the small-disturbance region. Eventually, the distinct 
vortex structures merge or break up into random, three-dimensional motions. 

Figure 1 illustrates the periodic vortex motion of a two-dimensional wall jet 
flow. Figure 1(a) shows a phase-averaged schlieren visualization obtained by 
lightly heating the nozzle flow. Figure 1(b) shows the phase-averaged velocity 
field, referenced to the convecting vortices, and Figure 1(c) depicts the corre- 
sponding vorticity field obtained with a central-differenced curl of the velocity 
field. The vorticity is normalized with respect to the maximum exit velocity and 
the nozzle width. The measurements were obtained from a single x-wire velocity 
probe by sampling the probe output at regular phase intervals as determined by a 
fixed pressure-transducer in the wall (ref. 17). These measurements were further 
processed to obtain estimates of the mean and fluctuating components of the 
dissipation field. 
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Figures 2, 3, and_4 show contour plots of the respective distributions of 
square mean velocity, o)_, mean square of fluctuating vorticity, u' , and total mean 
square vorticity, u + co' . Note that only the periodic component of the fluctuat- 
ing vorticity is resolved with the measurement scheme. The schlieren photograph of 
the flowfield in figure 1 reveal the existence of irregular, convecting vortex 
structures beyond the region indicated by the measurements. 


An overall view of the dissipative structure of the wall jet is gained by 
examining the variation of the dissipation integral with downstream distance. Here: 


where 


and 



(27) 


Distributions of and I2 are shown in figures 5 and 6, respectively. The value 
for !■] is nearly constant throughout the wall Jet until the end of the small- 
disturbance region, where it decreases by about 25?. The value of returns to 

about its initial level within the large-disturbance region. The overall vortical 
dissipation within the measuring region is approximately 1.1? of the kinetic energy 
supply from the nozzle. The total-acoustical-power output of the flow was estimated 
from reference 18 to be 0.02? of the kinetic energy rate. 


These measurements demonstrate a potential method of directly measuring dissi- 
pative terms in an unsteady flow. A direct extension of this method may be applied 
to experimentally determine a relationship between overall dissipation and variable 
parameters, such as forcing frequency, and to experimentally search for dissipation 
extrema. 


CONCLUSIONS 


An exact equation has been derived for the dissipation function of a homoge- 
neous, isotropic, newtonian fluid, with terms associated with fluid dilatation, wave 
radiation, and square vorticity. By using entropy extremum principles as given by 
Ziegler, simple flows such as the incompressible channel flow and the cylindrical 
Lamb vortex are identified as extremal configurations. The implications of the 
maximum dissipation conditions are found to be consistent with the principal notions 
of stability of parallel shear flows. The measurements of a typical periodic shear 
flow, the rectangular wall jet, show that direct measurement of the dissipation 
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terms is possible and that for this particular flow, the dissipation rate is nearly 
constant along the length of the Jet. 

Further experimental and theoretical analyses are required to determine the 
global validity of entropy extremum principles with regard to complex flows and to 
determine the applicability of time-averaged analyses. The limited observations 
obtained thus far suggest further applications to the analysis of turbulent, sound- 
generating, and shock-producing flows. 
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a. 




Figure 1.- Phase averaged velocity field and flow visualization, (a) Phase-averaged 
flow visualization; (b) velocity vectors (relative to convecting vortices; and 
(c) vorticity contours. Conditions: nozzle width, h = 0.508 cm; nozzle aspect 

ratio = 20; wall length, L = 3.81 cm; parabolic velocity profile at nozzle exit 
with maximum exit velocity, = 13.85 m/sec; tone frequency, f = 600 Hz. 
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Figure 2.- Contour plot of the square of the mean vorticity distribution, oi h /U^. 
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Figure 4.- Contour plot of the mean-total-square vorticity distribution, 

2 2 2 
(lo) h /Uq. 
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Figure 5.- Distribution of the first dissipation integral: I^(x/h) vs. x/h, where 



The nonzero vorticity region is enclosed by yi < y < y2* 
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Figure 6.- Distribution of the second dissipation integral: l2(x/h) vs. x/h, 
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